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Lagragian Interpolation lends itself nicely to machine 
calculation, but when the indices are equal the missing 
terms of the products require special treatment o The 
following approach decreases the number of multiplications 
and divisions required and practically eliminates the 
necessity of checking for index equality. 



2. METHOD 



If we defihe 

P (X) = (-^"^-i) (X-Xp)....(X-X, ),...(X-X t) (X-X )_ 

then Lagrange's interpolation formula Qlj becomes 

(1) y=^ ^k^W 



n 



U (X-X.), 



i=l 
ij=k 



k^l 



RECOMP TECHNICAL BULLETIN NO. 7 (Cont.) 



Assume X ^ X, for all k"= 1, 2, •.«, n, (if X=X, , y'^j,)» 
We may then write 



n_ -f|-(x-x.) 

^T '^k I'^l 1 or, 

^^^^ (X-X, ) P, (X. ) 

k«l ^ k k k 



y* 



It'"' 



J 



n 



y^ 



k=l (X-X^)P^ (1^) 



or. 



(2) 



Write 



n 



i»i 



n 



z 



y,. 



k=l (Xj^-X) Pj^ (x^) 



The denominator term (X,^-X) of (2) is similar to the terms of 
?^ (X^), so let 

(X^-X^) (X^-X2)...(X^-X^.^) (X^-X) (X^-X^^^)...(X^-X„_^) (X^-XJ. 

(Since X, is needed as a constant minuend in forming P, (X, ) 

it can be replaced in the table by X when it is picked up, 
and restored after P^ (X, ) is formed). 
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V/e now have 



n 



y 



i'=l 



(X-X.) 
1 



n 



Yy 



k«l 



< (\) 



Assuming that the products are set originally to -1 for 1 \ 

i 

and +1 for the ?f_ this method requires (n+l) or n +2n+l 

2 
multiplications and divisions, as opposed to 2n for the 

conventional method. It also simplifies the "housekeeping" 
by eliminating the tests for i=k. 

(It is, of course, nicest if there exists an instruction to 
exchange the contents of the arithmetic register with the 
contents of a word in memory), 
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